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The concept of the macroscopic wave function is a key for understanding macroscopic quantum 
phenomena. The existence of this object reflects a certain order, as is present in a Bose-Einstein 
condensate when a single-particle orbital is occupied by a macroscopic number of bosons. We extend 
these ideas to situations in which a condensate is acted on by an explicitly time-dependent force. 
While one might assume that such a force would necessarily degrade any pre-existing order, we 
demonstrate that macroscopic wave functions can persist even under strong forcing. Our definition 
of the time-dependent order parameter is based on a comparison of the evolution of Wparticle 
states on the one hand, and of states with — 1 particles on the other. Our simulations predict the 
possibility of an almost instantaneous dynamical destruction of a macroscopic wave function under 
currently accessible experimental conditions. 
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I. INTRODUCTION 

Superconductors, superfluids, and atomic Bose- 
Einstein condensates are described in terms of a macro¬ 
scopic wave function, a notion originally conceived in 
London’s theory of superfluidity [l|: Instead of consid¬ 
ering the Schrodinger wave function ..., r^r; t) of a 

Bose-condensed interacting fV-particle system, one works 
with an effective single-particle wave function $(rU) 
which obeys the nonlinear Gross-Pitaevskii equation 
[^. Experimental justification for the concept of the 
macroscopic wave function is provided by the observation 
of Josephson tunneling M between two superconductors 
coupled by a weak link [8|. Moreover, the occurrence of 
vortices, as observed in a series of landmark experiments 
with Bose-Einstein condensates [MU , is a direct conse¬ 
quence of the existence of a macroscopic wave function. 
Obviously, the reduction of the full A^-particle dynamics 
to that of a single-particle wave function requires that 
the system under consideration is highly ordered. This 
order is connected to the idea that $(r;t) represents a 
macroscopically occupied single-particle orbital, so that 
the terms “macroscopic wave function” and “order pa¬ 
rameter” often are used synonymously [l^ . 

But now new experimental developments are posing 
new theoretical challenges. There is an increasing ten¬ 
dency to subject Bose-Einstein condensates to strong 
time-dependent forcing, so as to “engineer” novel sys¬ 
tems which may not be accessible without such forcing. 
For instance, dynamic localization and quasienergy band 
engineering has been demonstrated with Bose-Einstein 
condensates in strongly shaken optical lattices EH, 
and coherent control over the superfluid-to-Mott insula¬ 
tor transition has been achieved Moreover, giant 

Bloch oscillations have been realized with condensates 
in tilted, ac-driven optical lattices [13) [ill- Still further 
experiments have demonstrated time-reversal symmetry 
breaking in shaken triangular lattices , and controlled 


photon-assisted tunneling (^ . [2l| . A particularly ambi¬ 
tious line of this research addresses the realization and 
usage of tunable artifical gauge fields [23 - 1^ or, phrased 
more generally, the exploitation of Bose-Einstein conden¬ 
sates in strongly forced optical lattices for quantum sim¬ 
ulation purposes [^ . 


These activities lead to an important question: To 
what extent is the underlying order degraded if one sub¬ 
jects a macroscopic wave function to strong forcing? It 
has been emphasized already quite early that the so¬ 
lution to the time-dependent Gross-Pitaevskii equation 
does not represent a condensate if it becomes chaotic [131 ■ 
The obvious conflict between dynamical chaos and the 
possible existence of an order parameter has inspired fur¬ 
ther studies especially on (5-kicked condensates, both the¬ 
oretical and experimental ones [Isl - I^ . But while the 
nonlinear Gross-Pitaevskii equation naturally can pro¬ 
duce chaotic solutions in the presence of external forcing, 
the actual A^-particle system still is described by a linear 
Schrodinger equation. Hence, while the iV-particle wave 
functions cannot become chaotic in the sense of nonlinear 
dynamics, there should nonetheless be a certain quality 
of the time-dependent A^-particle system which decides 
whether or not the solution to the Gross-Pitaevskii equa¬ 
tion actually qualifies as a macroscopic wave function, 
and there should be a measure which quantifies the de¬ 
gree of order remaining in a Bose-Einstein condensate 
under the action of an external force. In this letter we 
suggest an approach to these issues which does not in¬ 
volve the familiar partitioning of the field o pera tor into 
a condensate part and a noncondensate part |32 - l^ , but 
focuses on the evolution of neighboring states in Fock 
space. This may be seen as similar in spirit to the char¬ 
acterization of the degree of chaos in classical dynamical 
systems by probing the way initially close trajectories 
separate in time. 
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FIG. 1. Scheme for constructing the time-dependent 
macroscopic wave function. An initial A^-boson state 
|^(to))iv develops in time according to the time-evolution 
operator U{t,to), giving |^(t))Ar. If one acts with the field 
operator ^{r) on the initial state and normalizes, one obtains 
subsidiary [N — l)-particle states |3'(r|to))iv-i, which also 
propagate in time. A candidate function <l?(r; t) then is in¬ 
troduced by taking the matrix elements of the field operator 
with and |5'(r|t)}iv-i. On the other hand, propagat¬ 

ing first and annihilating thereafter yields |H(r|t))jv-i. If the 
absolute value of the projection n—i («'(r|t)|H(r|t))jv-i equals 
unity to good accuracy, 3>(r; t) actually is a macroscopic wave 
function which obeys the Gross-Pitaevskii equation. 


II. THE ORDER PARAMETER 

A key question is how the time-dependent macroscopic 
wave function, if it exists, is obtained from the full N- 
particle state. A guide to the answer can be inferred 
from the discussion given by Lifshitz and Pitaevskii [s^: 
Let |'I'(t))Ar be a time-dependent A-particle condensate 
state, and let |l£'(t))Ar_i be a “like” state of — 1 parti¬ 
cles; then the macroscopic wave function, normalized to 
unity, should be given by 

d>(r;t)= lim /^/N , (1) 

AT —>00 

where ^(r) is the bosonic field operator. But this leaves 
open the question how to quantify the “likeness” of 
|^(f))Af and |d'(t))Ar_i; if two such states are “like” at 
one particular moment to-, they might not remain so un¬ 
der the influence of time-dependent forcing. Moreover, it 
seems desirable to get rid of the limit of an infinite par¬ 
ticle number, and to study the emergence of a “macro¬ 
scopic” wave function already when N is relatively small. 
With this background, we proceed as summarized by 
Fig. [H We start from an initially given A^-boson state 
I^(^o))a/, which is not necessarily a pure condensate. Un¬ 
der the influence of some force which does not need to 
be specified at this point it develops in time into the 
A-particle state |'I'(t))Ar, as determined by the system’s 
time-evolution operator U{t,to)- In order to generate 


suitable {N — I)-particle states for taking the matrix el¬ 
ements suggested by Eq. we act with the bosonic 
annihilation operators V’(^) on the initial state, and nor¬ 
malize the results, obtaining 


|^'(r|to))Ar-i 


^{r)\'i'{to))N 


( 2 ) 


These subsidiary states likewise evolve in time under the 
action of the very same evolution operator, giving states 
|l['(r|t))Ar_i. We then define a function <l)(r;t) according 
to 


VN^{r]t) = Ar-i(^'(r|t)|^(r)|4'(t))Ar . (3) 

Observe the difference to the above tentative prescrip¬ 
tion We employ not just one single subsidiary 

(N — l)-particle state, but infinitely many; in princi¬ 
ple, there is one state |lE'(r|t)) 7 v-i associated with each 
r considered. Still, d>(r;t) as defined by Eq. ([3]) should 
qualify as a macroscopic wave function for sufficiently 
large N, and obey the Gross-Pitaevskii equation, pro¬ 
vided the above “likeness”-condition is satisfied. This 
means that the state trajectories evolving from the re¬ 
spective initial states |'I'(r|to))Ar-i and I'll (to)) at in Fock 
space should not diverge from each other too much, in a 
suitable sense. To bring this intuitive idea into a precise 
form, we also annihilate a boson from the time-evolved 
A^-particle state, thus producing 

^{r)\'i'{t))N 

t))Ar_i = ^- 

^ ^(r)t/(t,to)|^(to))jv 

ll^(r)|'I^(i))^|| 

Then the scalar products 

R{r; t) = Ar_i()i/(r|t)|S(r|t))Ar_i (5) 

^ N{'^{to)\^''{r)U''{t,to)^(r)U{t,to)\^{to))N 

||^(r)|vI/(to))^|| ||V^(r)|^k(t))^|l 

have a particular significance: If |i?(r;t)| = I, the can¬ 
didate $(r;t) provided by Eq. ([3]) is a true macroscopic 
wave function, obeying the Gross-Pitaevskii equation. In 
general, the magnitude \R{r;t)\, varying between 0 and 
1, provides the desired measure of the degree of order of 
the time-evolving A^-boson system. 

The justification for this statement stems from the ob¬ 
servation that the proper macroscopic wave function has 
to satisfy the requirement 

Af|$(r;t)p = N{'i'i.t)\ip^r)^{r)\'ii{t))N , (6) 

demanding that its absolute square, multiplied by the 
particle number N, yields the exact A-particle density 
of the system [s^. Introducing the projection operators 

Qt = |S(r|t))7v-i Ar-i(S(r|t)| , 


( 7 ) 
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we have the obvious identity 

= N{'i'{tMHr)Qt^{r)\'if{t))N- ( 8 ) 

Now, if this projector (I7|) were equal to the projector Pt 
defined by 


Pt = |«'(r|t))Ar_i Ar_i(^'(r|t)| , (9) 


which, in its turn, would be the case if |^(r|t))Ar-i dif¬ 
fered from |4'(r|t))Ar_i by not more than a phase factor, 
we could deduce 

= N{'i’{t)\^'^{r)Pti}{r)\'i>{t))N, ( 10 ) 

from which the desired identity dHl) follows immediately, 
keeping in mind the definition ([S]). 

This reasoning deserves still more scrutiny. Namely, 
if \E{r\t))N-i indeed differs from |T(r|t)) 7 v_i merely 
by a phase factor, then is proportional to 

|4'(r|t))Ar_i, wherefrom one is led to the relation 

^(r)|^'(t))Ar = \/iV$(r;t)|$(r|t))Ar_i . (11) 


This is reminiscent of what defines a condensate: As¬ 
suming that the N-particle state at some moment to 
corresponds to an N-fold occupied single-particle orbital 
(p{r,to) and thus is a pure condensate state of the form 


|'I'y,(to))Af = 




1 N 


d^rip{r;to) 


|vac) , (12) 


it obeys the equation 

i’{r)\'i'^{to))N = '/N(p{r-tQ)\'i>y,{to))N-i (13) 

at that moment to- Therefore, the condition m, which 
(if satisfied) guarantees that the time-dependent order 
parameter takes on its maximum value |it(r;t)| = 1, 
and thus makes sure that the candidate <i)(r;t) defined 
through Eq. © actually is a macroscopic wave func¬ 
tion, generalizes the familiar characterization of a pure 
condensate expressed by Eq. m so as to also involve 
time evolution, and reduces to it when the time t is 
close to to- Again adopting the dynamical-systems view¬ 
point, the projection ([5]) compares the trajectory of the 
given A^-particle state in Fock space to that of subsidiary, 
neighboring {N — l)-particle states. If it does not mat¬ 
ter whether one annihilates first and propagates then, or 
whether one propagates prior to annihilating, the flow in 
Fock space may be considered as (locally) stiff. Hence, we 
refer to the magnitude \R{r;t)\ as stiffness, with maxi¬ 
mum stiffness \R{r;t) \ = 1 expressing time-preserved co¬ 
herence in the sense of Eq. (El)- Note that the formal 
employment of subsidiary (TV — 1 (-particle states is nec¬ 
essary only to provide a reference for the evolution of the 


true TV-boson system: We do not violate particle num¬ 
ber conservation, and hence do not involve spontaneous 
symmetry breaking [3^ - 1^ . Moreover, an interesting ob¬ 
servation can be made here: If the decisive relation ED 
is satisfied, then the definition ([5D immediately yields 


R{r;t) 


d>(r;t) 
|$(r;t)| ’ 


(14) 


meaning that the phase of R{r;t) equals that of $(r;t). 
Read in the reverse direction, this implies that the phase 
of a macroscopic wave function contains information on 
the difference of the evolution of “neighboring” N- and 
{N — l)-particle states. This is well known in the equilib¬ 
rium case, when the phase of the solution to the Gross- 
Pitaevskii equation is determined by the chemical poten¬ 
tial, i.e., by the energy required to add one more particle 
to the system. The present considerations show that the 
phase retains a similar meaning even in case of nonequi¬ 
librium, in which a chemical potential does not exist. 


III. NUMERICAL SIMULATIONS 


In order to illustrate some consequences of the con¬ 
cepts developed above, we utilize the model of a bosonic 
Josephson junction [^, as described by the Hamil¬ 
tonian 


Hn = - 


^<2102 + +hK (^ a \ a\aiai + 03020202^ . 


(15) 

Here the hopping matrix element between the two sites 
labeled 1 and 2 is given by Ml/2, so that Ml is the 
single-particle tunneling splitting, while 2hK quantifies 
the repulsion energy of two particles occupying a com¬ 
mon site. The bosonic operator annihilates a particle 


at the jth site; aj is its adjoint creation operator. This 
system is subjected to a time-dependent bias with carrier 
frequency w and envelope fi/r(t), as specified by 


Hi{t) 


hn(t) sin(wt) — 02 ^ 2 ) 


the total Hamiltonian then reads 


(16) 


H{t)=Ho + Hiit) . (17) 

Even with constant amplitude /r(t) = /ii this model cap¬ 
tures nontrivial features of many-body dynamics [40l |41| ; 
it is one of the very rare systems which allows one to 
monitor the emergence of a macroscopic wave function 
numerically, but without further approximations on the 
iV-particle level. For all following simulations we select 
the ground state of the time-independent junc¬ 

tion ED with scaled interaction strength Nk/H = 2.0 
as the initial state. Note that this ground state is 
no A^-fold occupied single-particle state in the sense of 
Eq. (I12L because the relatively strong interparticle in¬ 
teraction leads to sizeable depletion [^. The required 
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FIG. 2. Degradation of the order parameter. Shown 
are the stiffness |i?i| (above) and the scaled population im¬ 
balance (1201) (below) for the driven bosonic Josephson junc¬ 
tion (II3 with N — 100 particles and scaled interaction 
strength Nk/Q. = 2.0, responding to a pulse with carrier fre¬ 
quency cu/fi = 1.6 and Gaussian envelope (11911 with width 
a/T = 10 and maximum driving strength /Xmax/fi = 0.51. 
The time scale is given by the cycle time T = 2'n/uo. The ini¬ 
tial state was the ground state of the undriven junction (nsj. 
Observe that the macroscopic wave function remains well pre¬ 
served until the middle of the pulse, after which the decrease 
of stiffness signals its degradation. 


subsidiary states then are given by 



N-l 




(18) 


for j = 1,2. Moreover, we fix the scaled carrier frequency 
uj/fl = 1.6, and consider a Gaussian envelope 


= Aimaxexp(-t^/2CT^) (19) 


with width a/T = 10, where the time scale is set by 
T = 2^!^}. In Fig. [5] we monitor the response of a sys¬ 
tem with N = 100 particles to a pulse with maximum 
driving strength = 0.51 by plotting the scaled 

population imbalance 

{J,){t)/N = jv(^'(t)|4ai - 4a2|^(f))Ar/(27V) (20) 

vs. time. We also show the stiffness |i?i(t)|; the cor¬ 
responding quantity |i? 2 (i)| obtained for the other site 
looks practically identical. Although N = 100 is not 
“macroscopically large”, one observes that |i?i| stays 
close to unity almost until the pulse’s middle, and then 
decreases in an oscillating manner. Thus, already in this 
situation there exists a good macroscopic wave function 
during the first half of the pulse, but it degrades signifi¬ 
cantly during the second half. 

Increasing the particle number to iV = 1000, while 
keeping Nk/VL and all other parameters constant, we ob¬ 
tain Fig. [31 This is a truly remarkable finding: Although 


FIG. 3. Preservation of the order parameter. As Fig. O 
but with N = 1000. Here the macroscopic wave function does 
hardly degrade during the entire pulse. Observe the scale of 
the insets’ ordinate! 


the iV-particle state undergoes violent changes when ad¬ 
justing itself to the driving force, the stiffness remains 
close to its theoretical maximum during the entire pulse, 
indicating that one can subject a macroscopic wave func¬ 
tion to strong forcing almost without reducing its order. 

A quite different scenario is depicted in Fig.|4l Here we 
have increased the driving amplitude to ^max/f2 = 0.55, 
and consider both N = 1000 (upper panel) and N = 
10000 (lower panel). While we observe excellent stiffness 
during the first half of the pulse, with 1 — |i?i(t)| appar¬ 
ently scaling with 1/A^, the macroscopic wave function 
is destroyed suddenly; this sudden destruction cannot be 
counteracted by an increase of N [i^. We utilize this 
example also to illustrate one more feature: As long as 
it exists, the macroscopic wave function should conform 
to the Gross-Pitaevskii equation. One may still solve 
that equation even beyond the point of destruction of 
the macroscopic wave function, but then the solution no 
longer captures the actual 7V-particle dynamics. This is 
verified by Fig. [SI where we superimpose the A^-particle 
imbalance (EOl) for N = 1000 to the prediction made 
by the Gross-Pitaevskii equation. As long as there is 
close-to-perfect stiffness, both curves are almost indis¬ 
tinguishable from each other, confirming the accuracy of 
the Gross-Pitaevskii approach under conditions of time- 
preserved coherence. But when the macroscopic wave 
function is destroyed the Gross-Pitaevskii dynamics be¬ 
come chaotic, losing their connection to the A^-particle 
level. 


IV. DISCUSSION 

The observations made in this work have both con¬ 
ceptual and experiment-oriented consequences. We have 
addressed the question why the solution to the time- 
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FIG. 4. Dynamically induced destruction of macro¬ 
scopic wave functions. As Fig. O but with higher driving 
amplitude /Tmax/fl = 0.55, and N = 1000 (upper panel) or 
N = 10000 (lower panel). The macroscopic wave function is 
not destroyed gradually, but quite suddenly; this destruction 
cannot be prevented by increasing the particle number. 


dependent nonlinear Gross-Pitaevskii equation can pro¬ 
vide a good description of a forced condensate only when 
it behaves in a regular, non-chaotic manner: That dis¬ 
tinction between order and chaos should have a coun¬ 
terpart already on the linear A^-particle level. As one 
possible characterization of this difference we suggest to 
monitor the time evolution of “neighboring” trajectories 
in Fock space of states consisting of N and iV — 1 parti¬ 
cles, respectively. With iV-particle states being orthogo¬ 
nal to states consisting of one particle less, the required 
measure of proximity of these states is provided by the 
projection of the former after annihilation of one parti¬ 
cle onto the latter. In this way, one can not only give 
a more definite meaning to the sketch by Lifshitz and 
Pitaevskii on how to construct the wave function of the 
condensate [1^, but one also obtains the desired indica¬ 
tor for the quality of this construction: The presence of 
a time-dependent macroscopic wave function necessarily 



ijT 

FIG. 5. Gross-Pitaevskii vs. A-particle dynamics. The 

A-particle population imbalance (I20II for /imax/fl = 0.55 and 
A = 1000, already recorded in the upper panel of Fig. [T] 
is compared to the prediction of the Gross-Pitaevskii equa¬ 
tion. As long as the stiffness is close to unity, there exists 
a macroscopic wave function which is perfectly described by 
the Gross-Pitaevskii equation, so that both curves almost co¬ 
incide. When the macroscopic wave function is destroyed the 
solution to the Gross-Pitaevskii equation becomes chaotic, 
and does no longer predict the A-particle dynamics correctly. 


requires that initially close state trajectories stay close to 
each other in the course of time. If this condition is satis¬ 
fied, the unmodified Gross-Pitaevskii equation provides 
an excellent description of the A-particle dynamics; if 
not, the macroscopic wave function is destroyed (d^ . 

Our matter-of-principle discussion is of little practi¬ 
cal help when it comes to computing the instability of 
a driven Bose-Einstein condensate in experimentally re¬ 
alistic situations, implying that knowledge of the ex¬ 
act A-particle state cannot be obtained. In such cases 
one requires other approaches, such as the second-order 
number-conserving self-consistent treatment developed 
by Gardiner and Morgan [s^, which has been applied 
to a toroidally trapped, d-kicked condensate by Billam et 
al. [Hill One then couples the solution of a general¬ 
ized Gross-Pitaevskii equation to modified Bogoliubov- 
de Gennes equations, assuming that the ratio of non¬ 
condensate to condensate particle numbers be a small 
parameter. This approach allows one to assess driven 
condensate dynamics with experimentally realistic parti¬ 
cle numbers jssl - f^ . 

Yet, even our idealized model calculations, which are 
not tied to any small parameter, do convey messages of 
practical importance. We have shown that a driving force 
does not necessarily destroy a macroscopic wave func¬ 
tion when it is applied smoothly, in the form of forcing 
pulses with a sufficiently slowly changing envelope. This 
particular manifestation of the quantum adiabatic prin¬ 
ciple signals green light for systematic quantum engineer- 
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ing with macroscopic wave functions. The identification 
of maximum stiffness, or of time-preserved coherence in 
the sense of Eq. dm), as the salient feature of a time- 
dependent macroscopic wave function may guide future 
investigations. Our simulations also illustrate an impor¬ 
tant fact: The initial 7V-particle state considered therein, 
which is the ground state of the model dTSl) , equals a pure 
condensate state only for vanishing interaction, that is, 
for Nk/VI = 0 [12, whereas we consider strong interpar¬ 
ticle interaction, Nn/fl = 2.0. Nonetheless, maximum 
stiffness can still be attained to an amazing degree of 
accuracy, as exemplified in Fig. [31 Finally, the observa¬ 
tion that substantial degradation of the underlying order 
parameter may not occur gradually in time, but rather 
can be connected to certain critical driving strengths, is 
open to experimental verification. Such experiments do 
not necessarily require a driven bosonic Josephson junc¬ 
tion, but can also be performed in other configurations. 


For instance, one could subject a Bose-Einstein conden¬ 
sate in a strongly anharmonic trap to a smooth forcing 
pulse, and perform a tiine-of-flight measurement of the 
condensate fraction after the pulse is over. If one repeats 
this measurement with successively stronger pulses, one 
should observe a sudden disappearance of the condensate 
peak at a certain critical maximum driving amplitude. 
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